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Abstract. We derive various approximations for the solutions of nonlinear 
hyperbolic systems with fastly oscillating initial data. We first provide er- 
ror estimates for the so-called slowly varying envelope, full dispersion, and 
Schrodinger approximations in a Wiener algebra; this functional framework 
allows us to give precise conditions on the validity of these models; we give in 
particular a rigorous proof of the "practical rule" which serves as a criterion for 
the use of the slowly varying envelope approximation (SVEA). We also discuss 
the extension of these models to short pulses and more generally to large spec- 
trum waves, such as chirped pulses. We then derive and justify rigorously a 
modified Schrodinger equation with improved frequency dispersion. Numerical 
computations are then presented, which confirm the theoretical predictions. 



1. Introduction 

1.1. General setting. The nonlinear Schrodinger equation has been derived as an 
asymptotic model for many physical problems involving the propagation of slowly 
modulated oscillating plane waves. Typical examples are water-waves |16[ [24' or 
ferromagnetism 20j and of course nonlinear optics (see e.g. [IHIEI]) for which it 
plays a central role. It is aimed at approximating the solution it of a nonlinear 
hyperbolic system with fast oscillating initial condition, say, 

f dtu + A{d)u + -Eu = eF{u), 

here e <C 1 is a small parameter corresponding to the wavelength of the oscillations 
in dimensionless variables (see Assumption [T] below for more precisions). The so- 
called Schrodinger approximation can be decomposed into two steps: 

-a- Slowly Varying Envelope Approximation (SVEA): one writes the solution 
u as the product of a fast oscillating wave train and an envelope U: 

(2) u{t, x) - U{t, a;)e' '""'^^""' + c.c, 

where (w(k),k) solves the dispersion relation (or, equivalently, belongs to 
the characteristic variety, see (fT3)) below), 
-b- The envelope U{t,x) is approximated by the solution of the nonlinear 
Schrodinger equation 

(3) dtU +{c^-V)U -e'-n{d,d)U ^eF{U), U\^^, = U\ 

where the group velocity Cg, the second order differential operator TZ{d, d) 
and the nonlinearity F can be explicitly given in terms of the data. 
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The Schrodinger approximation ©-Q has been rigorously justified [M l fTS t flT] 
for times of order 0{l/e) in the usual situation where the typical scale for the 
space variations of C/" is of order 1. This means that the number of oscillations (or 
optical cycles) in the laser pulse is of order 0{l/£). Recently, however, lasers with 
ultrashort pulses have been developed, for which the number of optical cycles is 
much smaller; for such pulses, the Schrodinger equation proves completely inaccu- 
rate and various authors proposed other ways of describing the asymptotics of ^ 
when e ^ 0. Some of these results are briefly recalled below; their common point 
is that they all abandon the SVEA, because the widely accepted "practical rule" 

(4) The SVEA Q is valid if |VC/°|oo < 7 
is enforced when the pulses get very small. 

Alterman and Ranch [P, f^j modeled short pulses by replacing the fast oscillating 
term in the initial condition by a fast decaying one; more precisely, they modified 
the initial condition for ([1]) as follows: 

(5) «|,^„ = [/°(xK^+c.c. - u\,^,,^U\x,^), 

with z) — > as z ^ 00, and the SVEA ^ is consequently replaced by 

u{t, x) - Uit, X, ^) , 

e 

with ?7(i, x, z) ^ as z ^ 00. The Schrodinger equation ([3]) is then replaced by 

(6) dtd.u + {c^-v)d,u + £]^n{d,d)u ^ed,F{u), c/u=o = t^°; 

this approximation (rigorously justified) uses the fact that the group velocity Cg 
does not depend on |k| and is therefore only valid in nondispersive media (E = 
in ([1])). Alterman and Ranch's approach has been generalized in [131 US] taking 
into account the particularities of the optical susceptibility of some cubic nonlinear 
media such as silica, and finally obtaining a quasilinear variant of ©J which is 
rigorously justified in the linear case. 

In order to model the propagation of ultrashort pulses in dispersive media, Barrailh 
and Lannes ^ chose another approach based on the functional tools developed in 
[T8j, which consist in replacing the initial condition for ([T]) as follows 

= C/°(a;)e*^ +C.C. = [/°(a;, 0, ^^), 

where the Fourier transform of the initial profile U^{x, T, Z) with respect to T and 
Z is an H'^ -valued measure of bounded variation. This general framework allows one 
to consider initial data of the form ([!]) -for which the bounded variation measure 
is obviously f°(a;)5(j^(k),k)~ and of the form ([5]). A generalization of ([6]) is then 
derived and rigorously justified; this equation however has the drawback of being 
quite complicated because the transport operator 3f + Cg • V must be replaced by 
a nonlocal operator modeling the fact that the group velocity Cg depends on the 
frequency in dispersive media. 

An important characteristic of ultrashort pulses that we did not mention so far 
is that their frequency spectrum is broad (while for usual wave packets as the initial 
condition of ([T]), and taking the Fourier transform with respect to the fast scale x/e, 
it is essentially contained in a 0(e) neighborhood of k; see for instance [4]). Since 
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the dispersion relation of the Schrodinger equation Q is a second order Taylor 
expansion of the exact dispersion relation of ([T]) at k, the error is quite important 
if frequencies far from k must be taken into account; in addition to the violation 
of the practical rule (U, this is another reason why the Schrodinger approximation 
breaks down for short pulses. This phenomenon is not specific to short pulses since 
it occurs for all pulses with large frequency spectrum (typical examples are chirped 
pulses). 

An alternative way to replace the NLS approximation for such pulses is therefore to 
focus on the dispersive properties of the asymptotic model. Instead of abandoning 
the SVEA (l2l)asin[Tll2[3[23l[13l[6], various authors [TTl[12] chose to make this 
approximation but kept the full dispersive properties of the original equations p]). 
thus avoiding the disastrous (for large spectrum pulses in dispersive media) second 
order Taylor expansion of the dispersion relation. Consequently, the approximations 
thus obtained might have a slightly smaller range of validity, but are undoubtly 
simpler and moreover provide a more precise approximation {0(e) versus o(l)). 
For instance, the equation derived in pjj, and which we call full dispersion model 
here, reads 

dtU+-(uj(k + eD)^uj.)U = eF{U), U{x) = f/°(x), 
e 

where (jj{-) parameterizes the graph of the relevant sheet of the characteristic variety 
(quite obviously ^ can be deduced from this equation by Taylor expanding uj{-) 
at k). This model is much simpler than the one derived in 6] and furnishes very 
satisfactory results; however it is still a nonlocal equation and its resolution requires 
spectral methods. One of the goals of this paper is to derive a new approximation, 
with, in practical, the same dispersive qualities as the full dispersion model, but 
keeping the same level of complexity as the usual Schrodinger equation dS]). 

Before describing with more details the results of this paper, let us introduce 
here two kinds of initial data for which we will often refer to throughout the 
article: 

• Short-pulses: the initial profile C/°(a;) in ([T|) is taken of the form 

(7) [/0(x)=/(^), 

with < /3 < 1 and / a smooth function; the case (3=1 corresponds to 
classical laser pulses, and short pulses to /? ^ 1 (see Fig. [Ij; the number 
of optical cycles for such a short pulse is thus 0{e/ fi); 

• Chirped pulses: the initial profile U^{x) in ([1]) is taken of the form 

(8) u\x) = f{x-xo)cos{-cos{^^)), 

P P 

with < /? < 1 and / a smooth function; the case /3 = 1 corresponds to 
classical laser pulses, and chirped pulses to /? <C 1 (see Fig. [2|). 

More precisely, we propose here to: 

(1) Provide a framework simpler than 6 but general enough to handle any kind 
of large spectrum pulses (such as ultrashort and chirped pulses). The idea 
here is to keep the SVEA and to work with envelopes which are in a Wiener 
algebra. The reason of this choice is that when /3 — > 0, the Wiener norm 
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Figure 1. Initial condition U°{x)e^^ +c.c. with [7° as in ([7]) and 
/ = e-^', xo = 5, e = 0.01, and for /3 = 1 and /3 = 0.1 




Figure 2. Initial condition U'^{x)e^^ " +c.c. with t/° as in ([8]) and 

/ = e^^ , xq = 5, e ~ 0.01, and for /3 = 1 and f3 = 0.1 



(which controls the i°°-norm) of the initial envelopes ([7]) remain bounded 
while any Sobolev norm controlling the i°°-norm grows to infinity; 

(2) Rigorously prove the "practical rule" (gl). We show that the SVEA Q 
makes sense if |VC/°|vi/ ^ where | • li^^ is the Wiener norm (see 
below); in the case of short-pulses ([7]), this condition is equivalent to e <C /3; 

(3) Establish precise error estimates for the full dispersion and Schrodinger 
models. This allows us to give precise estimates on the range of validity of 
these models; 

(4) Derive and rigorously justify a new family of Schrodinger equations with 
improved frequency dispersion and whose formulation is purely differential 
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(without nonlocal operator). The idea is to approximate the nonlocal op- 
erator of the full dispersion model by a suitable rational function, following 
an idea which proved very useful in water-waves theory (derivation of the 
BBM equation from the KdV equation [7j , derivation of Boussinesq models 
with improved frequency dispersion [HI [9]). These equations read 

{l-ieh-V - e^v • B\7)dtU + (cg • V)C/ 

- e^nid, d)U + ieV ■ (Vwi(k)b^)V[/ + e^CiV)U = F{U), 

where b e C*, B e Mdxdi^) and C : x C'' x ^ C is a trihnear 
mapping (taking b = 0, B — and C = gives therefore ([3])); 
(5) Present some numerical comparisons between the different asymptotic mod- 
els, for ultrashort and chirped pulses. 

1.2. Organization of the paper. Scction[2]is devoted to the derivation and jus- 
tification of the different asymptotic models. Basic assumptions, definitions, and 
tools are first given in ^2.11 and we then proceed to prove the envelope approxima- 
tion ^ in W2.21 in particular, precise error estimates are established which allow 
one to give a rigorous proof of the practical rule ([3]). In W2.'S[ we then derive and 
give precise error estimates in a Wiener algebra setting of the full dispersion model 
( ^2.3.ip and the Schrodinger model f ii2.3.2p . Finally the new Schrodinger equations 
dS]) are derived in ^2A3l 

In Section [3l we present numerical computations in order to make a comparison 
of the different asymptotic models. 

1.3. Notations. - We denote generically by Cst a constant whose value may 
change from one line to another; 

- We use the abbreviation c.c. for "complex conjugate", so that for all a e C, 
a + c.c. — a + a] 

- If a, b e then a • b denotes the scalar (not hermitian) product: a • b = 

- For all X G M'', we write x = {xi, . . . ,Xd), and dj stands for dj = dx ; we also 
write V = (^1, . . . , ddf and, for aU muhi-index a € N'', d" = d^^ . . . d^] 

- For all j = 1, . . . , d, we write Dj = Idj and D = (Di, ...,Dd)^ iV^; 

- The Fourier transform of a distribution u E S(M'') is denoted either by !Fu or u; 

- We use the classical notation for Fourier multipliers: 

-life L°-iR''), we simply write |/U = l/U-. 



2. Asymptotic results 

This section is devoted to the study of the asymptotic behavior when s goes to 
zero of the solution to the initial value problem (ivp) 

dtU + A{d)u H u = eT{u,u, u), 

I 

where A{d), E and T(-, •, •) are defined in the next section. 
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2.1. Basic assumptions and tools. We make the following assumption on the 
operators A{d), E and T(-, •, •) which appear in pO|) : 

Assumption 1. i. The system is symmetric hyperbolic in the sense that for 
some n>\: 

• One has A{d) — X]j=i ^j^j' '^^^ '"^s n x n symmetric, real-valued, 
matrices; 

• The n X n matrix E is real and skew-symmetric; 
ii. The mapping 

is linear with respect to u^, and . 

Example 1. A standard model for the propagation of a beam in a Kerr medium is 
the Maxwell- Lorentz system which can be written in dimensionless form as 

dtB + curl E^O, 
dtE - curl B + iQ = 0, 
5*Q-i(i?-P) = |PpP, 

where (E,B) is the electromagnetic field, P the polarization and Q = edtP ■ This 
system is of the form HIC^) and satisfies Assumption\^ with n — Ad and 



(11) 



A{d) = 



E = 



Vx \ 

-Vx 



/ 

(the entries in the above matrices are dxd matrices); denoting 
(j = 1,2,3J, the nonlinearity is given by 



V 






V 






'Idxd 





^ 

Idxd 

Idxd 
-Idxd / 

= {Bi,E^,Q^,P^) 



I 



T{u\u^,u^) 











Example 2. A simple toy model is the following Klein- Gordon system 



(12) dtu + 





V 



-V 




-V 




with u-.R^ xRi^^ Ci+'^ and v e M'*\{0}. 

Quite obviously, U^) is of the form ilU\) and satisfies Assumption{l\with n — l + d, 



Aid)^ 

and T{u^,u'^,u^) = {u^ ■ u^) ^ 





V 

-v^ 
V 



E = 







Under Assumption [l] the matrix ^(k) + -f — X]j=i ^j^j + T hermitian for 
all k e M'^, and thus diagonalizable, with real eigenvalues. We can therefore define 
the characteristic variety as 



(13) 



C:={(c.,k)e 



E^ 



, uj is an eigenvalue of A(\l) H }; 

i 
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introducing 

(14) £(w,k) :=-co/ + A(k) + -, 

I 

one can equivalently define C as the set of aii (w, k) e M x M'* sucli that k) is 
not invertible. 

It is classicai and not restrictive for our present concern to make the following 
assumption on C: 

Assumption 2. There exists m G N and m different smooth functions ujj e 
C°°(R''\{0}) (j = l,...,m), such that for all k e M.'^\{0}, the eigenvalues of 
yl(k) + M- are exactly cdj(k) (j — I, . . . ,m). 

Under this Assumption [21 one can write, for all k e M'*\{0}, 

„ m 

(15) A(k) + - = ^^,(k)^,(k), 

where 7rj(k) denote the eigenprojector associated to '^j(k) (in particular, tTj e 
C°°(R'*\{0};A^„(C)). 

We finally need a last assumption on the wave number k of the initial data of 
the ivp IjlOp in order to justify the asymptotic equations derived in this article. 

Assumption 3. One has k 7^ and, with co = ti>i(k); 

• One has (3w, 3k) ^ C; 

• With the notations of Assumption\^ 

3co > 0, Vj~2,...,m, inf Icj — Ci;,-(k)| > cq. 

kGB<i\{0} ' ' 

Remark 1. The first part of the assumption excludes resonances with the third har- 
monic. The results presented here could easily be extended to cover such a situation, 
but it is not restrictive at all to make this assumption. 

Example 3. For the Maxwell equations one can check after some computa- 

tions that Assumptions^ is satisfied with m = 7 and 

cu,{k) ^ i(V2(l + |k|) + |kP + V2(l-|k|) + |kP), 
W2(k) = V2, 

c.3(k) = i(v/2(l + |k|) + |kp- V2(l-|k|) + |kp), 
W4(k) = 0, 

and Lo^ = — W3, loq — —L02, ^7 = —u>i (so that lui > uj2 > ■ ■ ■ > oj^ on M'^\{0}J; one 
can therefore take cq — wi(k) — a/2 > 0. 

Example 4. For the Klein-Gordon system one readily checks that Assumption 
is satisfied with m = 2, u;i(k) = y^kp + |vp and UJ2 — —uji. One can then 
remark that As.sumption\3[ also holds for allk.^0 and cq = LOi{h) + |v|. 
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We finally end this section with some results on the Wiener algebras. First 
recall that W''{R'^- C") (fc, neN) and W{Rj^ x T9; C") (which will be denoted by 
^'^ X T;C") in the sequel) are defined as 



(16) W\R'^; C") := {/ G 5'(R'^)",Va e N'^, \a\ < k, \d''f\w < 00}, 
with \f\w ■■= |/|li(r<';C"), and 

WiR" X T; C") := {/ = E l/k(«<^xT) := E < 

(when fc = 0, we write W{R'^; C") instead of VF0(R''; C")). 

The classical properties of the Wiener algebras used in this article are recalled in 
the following proposition. 

Proposition 1. i. The space W''{R'^;C), k e N, (resp. VF(R'^ x T;C); is an 
algebra in the sense that the mapping {f,g) 1— > fg is continuous from W''{R'^;C)'^ 
into W''{'^'^;C) (resp. W{R'^ x T;C)^ into W{R'^ x T;C);. 

ii. If M is a skew- symmetric, real valued n x n matrix, then exp(—iM) is unitary 
onW''iR'^;C") (resp. iy'=(R^ x T; C")/ 

iii. /// € W{R'^;C) then for all /3 > 0, one has /(^) e W(R'']C) and \ fi^)\w = 
\f\w- 

Remark 2. As said in the introduction, the third point of Proposition{J\is the main 
motivation to work with Wiener algebra rather than Sobolev spaces, because the 
W{R'^;€.''^)-norm of initial conditions of the form ^ remains bounded (constant) 
while its H^-norm is of size 0{(3'^/'^^^) as /? — )■ 0. This framework is somehow a 
simplified version of the functional setting of [S] (where the SVEA is not made) and 
has also proved useful in the study of wave-packets [5] or localized solutions [22j . 

2.2. The envelope approximation. This section is devoted to the justification 
of the envelope approximation which states that the exact solution to pI7)) can be 
described at leading order by u^pp defined as 

(17) ulpp{t,x)^U{t,x)e'^^ +C.C., 

where the envelope U solves the envelope equation 

dtU + ^£(w, k + eD)U = eT{U) 
U\,^„{x) = U^{x), 



(18) 

with £(•,•) given by (I14p and T is defined as 



T{U) = T{U, U, U) + T{U, U, U) + T{U, U, U). 
The interest of the envelope equation (fT5|) with respect to the original ivp H 



is that the fast oscillating scale has been removed from the initial data. The main 
result of this section is the following theorem. 

Theorem 1. Let AssumptionsUl [H and[3l be satisfied, and let a,b G W^{R'^; C"), 
and 

~ 7ri(k)a + eb (Polarization condition). 

i. There exists a time tq > such that for all < e < I, there is a unique solution 
U e C([0, f);W{R'^] C")) to m 
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ii. For all < T < tq, there exists sq such that for all < e < eo, there is a unique 
solution e C([0,r/e] x M'')" to i l70|) . and one has 

\<x - <pplL~([0,r/e]xR^)" < eC{T, \U''\w){l + \b\w + |V[/°|w), 

where M^^p is as defined in |J7| ). 

Remark 3. Assuming that \U^\w and \b\w o-re 0(1) quantities -which is of course 
the case for short pulses-, one deduces the following result 

The SVEA is valid if |V?7°|w < ^, 

and Theorem]^ thus provides a rigorous basis for the "practical rule" When 
working with short pulses with initial condition Qj, it is easy to check that this 
condition reads simply e <ti f3. 

Remark 4. Working in the more classical framework of Sobolev spaces, one could 
establish an error estimate similar to the one given by the theorem, but with H'^- 
norms (s > d/2) instead of W -norms in the rhs of the estimates. For short pulses 
with initial data the control would therefore be of the form eC (r, ^^s-d/i ) ps+(i-d)/2 
which is obviously useless when /3 — > 0. 

Proof. Let us prove the following lemma, which implies the first point of the theo- 
rem. 

Lemma 1. Let [/" g VF(M'';C"), and assume that Assumption[J\is satisfied. 
There exists tq > such that for all < s < 1 one has a unique solution U G 
C([0, f); WiW^; C")) to (Qg). For all < t < tq, one has 

(a) sup \U{t)\w <C{t,\U°\w); 

0<t<T/e 

if moreover € W^{R'^;C^) then one also has 

(b) sup |Vt/(i)k <C(T,|[/"|H.)|Vf/°|w'. 

0<t<T/e 

Proof. Uniqueness is obvious, and to prove existence, we use a classical iterative 
method: let Uq = U° and, for aU n G N, 

Un+i{t) = 5(t);7° + e [ S{t- t')T{Un)dt', 

with S{t) = exp(-f/:(w,k + ££>)). 

Since C{uj_,'k + e£^) is real and skew-symmetric for all ^ e M'', we can use Proposition 
[T] and the trilinearity of T to get 

sup\Un+i\w < \U°\w +£Cst t{swp\U„\wy, 

[OA] [0,t] 

and 

SUp|[/„+l - Un\w < eCst tC{ sup \Un\w, sup \Un~l\w) SUp | L/„ - Un-l\w, 
[0,t] [0,t] [0,4] [0,4] 

and it follows easily that the sequence (C/„)„ converges in C([0, |-]; M^"(M'^; C")) 
(for some r > 0) to a solution U of 

(19) U{t) ^ S{t)U" + e f S{t-t')T{U)dt'; 

Jo 
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the solution can then be extended to a maximal time interval [0, Isi^isi) and one 
can show with classical arguments that tq := info<e<i Tniax(£) > 0. The estimate 
(a) of the lemma also follows easily from (fT9|) and a Gronwall-type lemma. 
Differentiating (fT9|) with respect to Xj [j = and using Proposition [T] to 

control the VF(R'^;C") -norm, one gets 

\d,U {t)\w <\djU°\w + eCsi j \U{t')\lr\djU{t')\wdt' , 

Jo 

and the estimate (b) follows from Gronwall's lemma and the estimate (a). □ 

Before going further in the proof of the theorem, let us introduce some notation. 
We decompose the solution U of provided by Lemma [T] as 

[/ = C/i H h [/„ , with Uj = TTj (k + eD)U, 

and we also write Uu = U2 + ■ ■ ■ + Um- 

The first step of the proof of the second part of the theorem consists in controlling 
dtUi uniformly in e -which is much better than the 0{l/e) estimate on dtU one 
can deduce directly from the equation (fTS]) . 

Lemma 2. If Assumptions]^ and\^ are satisfied, then for all < t < tq, 
sup \dtUi\w < C(t, |C/°|m/)(1 + |VC/°|m/). 

0<t<T/e 

Proof. Multiplying the envelope equation p8|) by 7ri(k + eD), one gets 
(20) dtUi + ^{uJi{k + eD)^ui)Ui = STri(k + eD)T{U). 

Recalling that u = uji(k.), a first order Taylor expansion shows that 

\{uji{k + eD) - ui)Ui\^ < e\VLJi\Lo^\VUi\w 

(note that one infers |Va;i|Loo < 00 from the observation that for all k 7^ 0, 
9jtt'i(k)7ri(k) — 7ri(k)74j7ri(k)). It follows therefore from ([^0]) . the trilinearity of 
T, and Proposition [1] that 

\dtUi\w < |Vc^i|ioo|V?7i|M/ + Cste|t/|^, 

and the result follows from Lemma [TJ □ 

We now prove that the components Uj {j > 2) remain of size 0{e) if this is 
initially the case. 

Lemma 3. // Assumptions\Tl [H anc![3 are satisfied, and if C/" = 7ri(k)a + eb, then 
one has, for all < t < tq, 

sup \Uii{t)\w < eC{T, |(7°|w)(l + \b\w + \yU°\w). 

te[0,T/e] 

Proof Multiplying ^ by TTj{k + eD) {j > 2) gives 

dtUj + ^ {toj (k + eD) - ui) Uj = en, (k + eD)T{U) 

= enj (k + eD)T{Ui) + en, (k + eD) {T{U) - T{Ui)) . 
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With Sj (t) = exp ( — (k + eD) ~ uj_)), one gets therefore 

Uj{t) = Sj{t)U^ + e f Sj{t-t')TTj{k + eD)T{Ui)dt' 
Jo 

(21) +e [ Sj{t-t')n,{^ + eD){T{U)^T{Ui))dt'. 



We now bound the ly-norm of the three terms of the r.h.s. of pT|) : 

• Estimate of Sj{t)U^. Since Sj{t) is unitary on VF(R'*;C"), one has 



w- 



Since moreover one can write 

7rj (k + eD)U° = (7rj(k + eD) - 7rj(k))C/° + 7rj (k)C/", 
it follows from the orthogonality of the projectors tTj (j = 1, . . . , m) that 

TTj (k + ei:>)C/° = (ttj (k + eD) - tt^ (k));7° + ettj (k)6. 

Since the derivatives of 7rj(-) are in general not bounded near the origin, 
we cannot control the first term of the rhs by a Taylor expansion and we 
thus write 

^,(k + £i^)C/° - (^,(k + £Z?)-^,(k))l{,|,,|<|k|/2}C/° 

+ (^j(k + el?) - ^,(k))l{,p|>|k|/2}C/° + £^j(k)6, 

where l{e|5|<|k|/2} = 1 if £|^| < |k|/2 and otherwise. 
Using the fact that 7rj(-) is C°° on the ball of center k and radius |k|/2, we 
can bound the first term of the rhs in Vl^(R'', C")-norm by eCst |V[/°|vi/; 
one can also check that a similar estimate holds for the second term of the 
Ihs since one has 1 < pfj-lCl for ^I'CI ^ lk|/2- We can thus conclude that 

(22) \Sj{t)U'^\w < eCst {\b\w + |V{/"|w)- 

• Estimate of A := e J* Sj{t - t')Trj{k + eD)T{Ui)dt'. Taking the Fourier 
transform of this term and integrating by parts yields 

/■* t -t' — - — 
e / exp(-i {u;j{k + e£_) - ui))nj{k + e^)T{Ui)dt' 

Jo £ 

• /•* £exp(-i^^(^,-(k + eO-^)) .An-fTTTw/ 
= —le / ^ — -T:j(k + e^)OtT[Ui)dt 



ujj (k + e^) — w 



+ie 



eexp ( — i^—^{ujj{k + e^) — w)) 



ujj (k + e^) — w 
One deduces therefore, using Assumption [31 that 



7r,(k + eOT(C/i) 



e 



2 



|3 



sup 1^4(^)1^^ < Cst — sup \Ui\y^ sup |(9tC/i|Ty + Cst — sup |C/l|^, 

te[0,T/e] Co [0,T/e] [0,T/e] Co [0,T/e] 

SO that, owing to Lemmas [T] and [H 
(23) sup \A{t)\w <eC{T,\U''\w){l + \W'\w). 

te[0,r/e] 
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• Estimate of S := £ Sj{t- t')TTj{k + eD){T{U) -T{Ui))dt' . First remark 
that owing to the trihnearity of T, one has for all t e [0, r/e], 

\T{U)it)-T{Ui)it)\^<Cst sup \U\^w\Uii{t)\w; 

[O.r/e] 

using Lemma [U we obtain therefore 

(24) sup \B{t)\w <eC{T,\U°\w) f \Uii{t')\wdt' . 

te[0,T/e] Jo 

It is now a direct consequence of (l211) j {j = 2,...m) and ((22l) - ((24l) that for all 
te[0,T/e], 

\Un{t)\w < e{\b\w + C{t, |(7°|w)(1 + |Vt/°|w) + £C(t, |C7°|h.) / \Uii{t')\wdt' , 

Jo 

and the result follows therefore from Gronwall's lemma. □ 

We are now set to conclude the proof of the theorem. We look for an exact 
solution Mgj. to pop under the form 

. , . , k • X — wt 

with Uex e W^(M'^ X T; C") itself of the form 

Uex{t,x,e) = Uapp{t,x,e) + ev{t,x,e), 

with Uappit, X, 9) = U{t, a;)e*® + c.c. and V bounded in W{W^ x T; C"). With U as 
given by Lemma [l] the equation that V must solve is 

dtV + -C{uj_Dg,kDe + eD)V = T(U. U, U)e'^^ + c.c. 
e 

~\~{T{Uq^PP -\- Uapp ~^ 1 Uapp ~t~ ^ {Uapp ^ Uapp^ C^app)) ■ 

Owing to the first part of Assumption [31 we can look for V under the form 

V{t,x,0) = Voit,x,e) +eVi{t,x)e'^^ +C.C., 
with Vi = —iC{3uj_, 3k)^^r([/i, Ui, Ui); the resulting equation on Vq is 

(25) dtVo + -C{ijDe,kDe + eD)Vo =^h+h + h, 

£ 

with 

h = {T{U,U,U)-T{UuU,,Ui))^^' + c.c. 
h = -e{dt + A{d))Vie^''' + c.c. 

I3 = {T {Uapp + SV, Uapp + SV, Uapp + eV) — T {Uapp, Uapp, Uapp)). 

Let us now bound Ij {j = 1, 2, 3) in W{R'^ x T; C") and for all t e [0, r/e]: 

• From Lemmas [1] and [3l one gets 

(26) |/i(OI WxT) < eC{T, |!7°|w)(l + \b\w + |V!7°|w). 

• From the definition of Vi and Lemmas [l][2l one has directly 

(27) |/2(i)|w(E-xT) < eC{T, \U"\w){l + m"\w). 

• From the trihnearity of T and Lemma [TJ one gets 

(28) 1/3(01 WxT) < £C{t, \U°\w){l + \Vo{t)\w + emt)\l, + e^\Vo{t)\l,). 
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By Proposition [U the semigroup S{t) — cxp(—i^C{uj_D0,'kDe + eD)) is unitary on 
X T), so that the estimates ([26|) - ([28l) ahow one to conclude to the existence 
of a solution Vb e C{[{),t / e\;W {W^ x T)") to ^ using a fixed point formulation 
similar to the one used in the proof of Lemma [1] After a Gronwall argument, one 
also gets 

(29) sup |yo(i)k<C(T,|C/0|w)(l + |fok + |VC/°|H'). 

0<t<T/e 

Remark 5. In fact, one finds that Vq exists a priori on a time interval [0,r'/e], 
with t' < T. However, by a classical procedure of continuous induction (such as 
in the proof of Theorem 3 in |19j for instance), one can get t' = t, provided that 
< e < £o with £o small enough. 

It follows from the above that 

sup \U,Ai) - C/appWlH/(R''xT) < eCir, |C/0|h')(1 + \h\w + |Vi7°|iy), 

te[0,r/e] 

and the theorem follows therefore from the observation that 

l^'L - "Sppl-L~([0,T/e]xR'i) < sup |t/ea;(i) - C/app(Oliy(R''xT)- 

te[o,T/£] 

□ 

2.3. Approximations by scalar equations. 

2.3.1. The full dispersion model. The full dispersion model consists in approximat- 
ing the exact solution to (fTO|l by u^app i defined as 

(30) <pp,i(*'^) = U^i){t,x)e'^^ +C.C., 
where C/(i) solves the full dispersion scalar equation 

r 9tC/(i) + f(^i(k + £i^)-c.)C/(i)-e7ri(k)T((7(i)) 
^ ' I U^^^u^,ix)=U'{x) 

and with wi(-) as in Assumption [T] 

The following corollary shows that the full dispersion scalar equation yields an 
approximation of same precision than the envelope equation for times t e [0, t/s]. 

Corollary 1 (Full dispersion model). Under the assumptions of Theorem\^ and 
for all < s < So (sq > .small enough), there exists a unique solution [/(x) G 
C([0,ro/e);W^(R'';C")) to fHP. 
For all < T < Tq, one also has 

\<x - <ppalL-([0,T/e]xR^) < £C{t, \U°\w){l + \b\w + |Vf/°|w), 
where w^^p ^ is as defined in \30\) . 

Remark 6. The quantity [/(i) remains C"-valued, but we call iSl]) a scalar ap- 
proximation because the operator | (t^^i (k + e_D) — tj) is scalar, which is not the case 
of ^>C(ti;, k + eD) in the envelope equation I118\) . The interest of the FD model is 
that ^{u!i{'k + eD) — uj_)u remains bounded for spectrally localized functions u, while 
|'C(<±^, k-l-eZ))^ is of order 0{l/e). The fast oscillations of the nonpolarized modes 
must therefore be taken into account with the envelope approximation, and the dis- 
cretization step must therefore be much smaller in numerical computations than for 
the FD model. 
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Remark 7. i. Performing the same analysis as in Remark\^ one can check that 
the "practical rule" also applies for the FD model. 

ii. The FD model has been derived and studied in '12] (it is called "intermedi- 
ate model" in that reference) for the study of chirped pulses with initial data (0)- 
However the error estimate provided in [12 is of the form £C(t, jt/^l^fs), with s 
large, which, as explained in Remark^is not small for short pulses (or more gen- 
erally large spectrum -including chirped- pulses). Note that a variant of the FD 
model with the orthoganal projector tti (k) replaced by tti (k + eD) in front of the 
nonlinearity is also studied in (12j and gives very good results. 

Proof. We omit the existence/uniqueness part of the corollary, since it is obtained 
with the same tools as for Theorem[T](in particular, taking a smaller eg if necessary, 
the existence time of the envelope equation is larger than the existence time for (j3ip 
and we can thus take the same tq as in Theorem [Ij, and we thus focus on the error 
estimate. 

Denoting as in the proof of Theorem [T] C/i = 7ri(k + eD)U, where U is the solution 
of the envelope equation, one gets from Lemma [3] that 

sup Pit) - Ui{t)\w < eC{T, |C/°k)(l + \b\w + |VC/°|w), 

te[0,T/e] 

so that it suffices to control \Ui{t) — to prove the corollary. Applying 

Tri{k-^eD) to one gets 

dtUi + ^{LOi{k + sD)^ui)Ui ^ eT:i(k + eD)T{U), 
so that the difference F = C/i — C/(i) solves 

f^r.. / dtV+^iuJi{k + eD)-iu)V = e7:i{k + eD)T{U)^e7ri{-k)T{U^r)) 

Remark now that 

TTi (k + eD)T{U)- TTi (k)r ([/(I) ) ^ {Tii{k + eD) - tti (k)) T{U) 

(33) +7ri(k)(r([/) - T{U,)) + 7ri(k)(r(t/i) - r({/(i))). 

Since |(7ri(k + ei:*) - T:i{k))r{U)\w < eCst \VT{U)\w (see the proof of ((22|)). 
one can use Lemma [T] to bound the first component of the rhs of ([33|) from above 
by eC(T, |J7°|vy)|V?7°|w. The second component of (|33|) can be estimated exactly 
as the term Ii in ((26|) . while the last one is bounded from above in W{M.'^;C") 
by C(|C/i|vk,|^^(i)|h')|^|h'- Since moreover \Ui\w is controlled by Lemma [T] and 
that a similar estimates also holds obviously for |J7(i)|w, one deduces that for all 
< < < r/e, 

\T{U{t))-T{U^,){t))\w <£C{T,\U''\w){l + \b\w + \^U''\w) + C{T,\U''\w)\V{t)\w. 

This inequality, together with an energy estimate on (|32p and a Gronwall argument 
shows that 

sup \V{t)\w < eC{T, |C/°|w)(l + \b\w + NU"\w), 

te[0,T/e] 

where we also used the estimate \TTi{k + eD)U° -U"\w < eCst (1 + + |VC/°|w) 
(which is proved with the same arguments as (|22p ). □ 
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2.3.2. The nonlinear Schrodinger equation. In the Schrodinger approximation, the 
exact solution to ()10p is approximated by w^pp 2 defined as 



(34) ulpp^it^x) ^U^2){t,x)e' - +C.C., 

where C/(2) solves the nonlinear Schrodinger equation 



(35) J 9*C/(2) + (Vc^i(k)-V)C/(2)-e^(V-i7k(a;i)V)[/(2)=e7ri(k)T(;7(2)) 

I f/(2)U„(x) = {7"(a;), 

and where H-kitoi) stands for the Hessian of ti;i(-) at k. One then has 

Corollary 2 (Schrodinger approximation). Under the assumptions of Theorem\^ 
and for all < e < Sq (eq > small enough), there exists a unique solution 
[/(2) € C([0,To/£);M^(R'^;C")) to 

If moreover E W^iW^; C") then for all < t < tq, one also has 

|wL-<pp,2lL~([0,r/e]xR<i) < >^C{t, \U°\w)il + l^U^'lw + \b\w + \cSchrod\oo\U"\w3) , 

where u'^pp 2 *s as defined in ^34^ and 

_ {uj.jk + ep - + £Vo7i(k) -i + em- gk(cc>i)e)) 

Remark 8. i. A third order Taylor expansion o/aji(k + e^) at ^ — shows that 
\^Schrod\oo is finite and can be bounded from above independently from e. 

ii. The component \cschrod\oo\U'^\w'^ of the error estimate does not appear for 
the full dispersion model. It is due to the approximation of the nonlocal operator 
7(wi(k + eD) - w) (Ihs of fSJ]) ) by the differential operator Vwi(k) • VC/(2) - e|V • 
_ffkV?7(2) (Ihs of !i35\} }. This error term is thus a linear effect. 

iii. This additional term is responsible for the bad behavior of the Schrodinger equa- 
tion to model short pulses ( and more generally large spectrum waves such as chirped 
pulses). For instance, for initial data like the precision of the Schrodinger 
approximation is of order 0(e(i + llsahrodU^^-^ when j3 <^ 1. In order for the 
Schrodinger approximation to keep the same order of precision as the full dis- 
persion model, one needs therefore to have ^'^^''^'3°'^^°° < jj, which requires that 

\(-Sch7-od\oo- This condition is far much restrictive than the practical rule 

/3 > £. 

iv. The Schrodinger approximation has been rigorously justified for systems of the 
form J JO)) ( and also for some quasilinear generalizations ) by Donnat, Joly, Metivier 
and Rauch \1A\ I15j and Lannes I17j but in the Sobolev framework which we saw is 
not adapted for the .study of large spectrum pulses. 

Proof. As in the proof of Corollary [l] we focus on the error estimate and omit the 
existence/uniqueness part of the proof. 

The difference = C/(i) — C/(2) of the solution of the full dispersion and Schrodinger 

equations solves the initial value problem 

(36) 

dtV + l{uJi{k + eD) - ui)V = e7ri(k)(r(C/(i)) - T(C/(2))) - e^TZ2{D)U^2) 
V\,^„ix) = 0. 

where for aU C G K'', 

7^2(C) = ^(zc^i(k + £0-i^-«£Vc^i(k)-C-e'^e-^k(^i)C)- 
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Remark now that one has for all < t < r/e, 

hrod\oo 

with Cschrodi') sls In the statement of the corollary; differentiating the Schrodinger 
equation l|55|) and estimating the W^-norm of the solution, one also gets easily 

sup |C/(2)(i)k3 < C(t, |C/V)(1 + \U°\w-^). 

telO,T/e] 

Since the first term of the r.h.s. of ([36| can be bounded as in (|26p. one gets from 
Gronwall's lemma applied to p6|) that 



sup \V{t)\w <CiT,\U°\w)il + \csc hrod\cfo\U Ixys) 

te[0,T/e] 

which, together with Corollary[l] yields the result. □ 

2.3.3. The nonlinear Schrodinger equation with improved dispersion relation. As 
said in the introduction, we propose in this paper new approximations based on a 
family of modified Schrodinger equations, whose dispersive properties are closer to 
the exact model. Such an approximation u'^pp 3 is defined as 

(37) <pp.,3it,x) = f/(3)(t,a:)e'^^ +c.c., 

where C/(3) solves the nonlinear Schrodinger equation with improved dispersion re- 
lation: 

[l-ieh-V - e^V • BV)dtU(^3) 

(38) ; +(Vc^i(k)-V-zeV-(ii7kK)+Vc^i(k)b^)V + £2c(V))C/(3) 
= e^i(k)r({/(3)) 

where b e C*, B G Mdy.d{^) and C : x x ^ C is a trilinear mapping. 
We assume moreover that 

(39) B is symmetric positive, b e Range{B), and 4 — b • (B^^h) > 

(note that even though B^^h is not unique when B is not definite, the scalar 
b • {B~^h) is uniquely defined). One then has the following result: 

Corollary 3 (Improved Schrodinger approximation). Under the assumptions of 
Theorem]^ and for all < e < (eq > small enough), there exists a unique 
solution U(3) e C([0,ro/£);H^(K'^;C")) to 

If moreover [/" G W'^{R'^; C") then for all < t < tq, one also has 

Wlx-Kpp,3\L°^{[0,T/e]xM'i) < sCiT, |f/° |h/) (1 + 1 VC/" | VF + 1&| 1^ + 1 Cimproued |oo | | H^a ) 7 

where u^^pp 3 is as defined in Jg?] ) and 

_ ^llk + eC) - + g i+^b-g+e^g-Bg ) 

timproved{C,} ^^3 ( _|_ |^ |3 ) 

Remark 9. i. As for the Schrodinger equation, one can check by a simple Taylor 
expansion that \cimproved\oo is finite and uniformly bounded with respect to e. 

ii. Taking b = 0, -B = and C = (this choice satisfies I139\) ). one recovers the 
usual Schrodinger equation i35\) . 

iii. The interest of i38\) with respect to i35\) is that one can choose h, B and 
C such that ^improved <C Cschrod, thus improving considerably the accuracy of the 
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Figure 3. 



for e — 0.01 with the coefficients (l40l 



approximation. In the one dimensional case d — I, it is possible to choose h, B and 
C in such a way that the dispersion relation for i38\) is the [3,2]-Pade expansion of 
the dispersion relation of \31\) . For the case of the Klein-Gordon system this 
leads to 

(40) b = ^, B.^^±i^, c.M^^:!±i^); 

^ ' v2+k2' 4(v2+k2)2 4(v2+k2)5/2' 

we illustrate in Figure [3 how much one gains by working with h3^) instead of k35fl 
for the Klein- Gordon system \l'/3fl with v = k = 1. 



iv. The same analysis as is Remark\^iii shows that the approximation provided 
by i38\) is of the same order as the envelope approximation if 0^ > \ Cr 
Since la, 



'^improved \ oo ■ 

vedloo ^ |cschrod|oo7 this Condition is much weaker than the correspond- 



ing one for the usual Schrodinger model. In some particular cases, this condition 
can even be weaker than the "practical rule" e ^ f3. 

Proof. Choosing G -~-^B^^h, one can check that 

1 + b • c + e • i?e = 1 - 7b ■ (B-'h) + (e - ^o) • B{( - eo), 



so that it follows from assumption ([55)) that 1 + b • ^ + ^ • > (uniformly with 
respect to ^ e M'*). The operator 1 — eih — e^V ■ BV is therefore invertible, and its 
inverse is the Fourier multiplier (1 + eb • _D + e^/) . BD)~^. The equation (|38p can 
therefore be rewritten as 



. VcJi(k) -D + eD- (ii?k(t^i) + \7uji{k)h'^)D - e'^C{D) 

(l + eh-D + e^D-BD) 
{1+eh-D + e^D- BD)-^tii{^)T{Ui^^)). 



(3) 



Since {l-\-eh- D -\-e'^D ■ BD)-'^ is regularizing (of order -2) and acts on H/(R'';C") 
uniformly with respect to e > 0, the proof of the result follows exactly the same 
lines as the proof of Corollary [3] and we thus omit it. □ 
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3. Numerical simulations 

3.1. The equations. This section is devoted to the comparison of solutions of the 
different asymptotic equations derived in Section [2] with the solutions of the full 
system (|10p. In particular, we want to check numerically the results proved in 
Corollaries [2 El and [3 

We consider both short pulses ( §3.3|) and chirped pulses ( §3.4|) : for the numerical 
computations, we use the toy model (I12p with dimension d — I and v = 1 (this is 
also the model used in jl2j for the study of chirped pulses). Writing u — {f,g)^, 
this model reads 

f dj + d^g-^^~e{\f\' + \g\')g, 

(41) < f 

[ dtg + d,f + ^=ei\f\' + \g\')f. 

The initial conditions are taken of the form 

(42) /|,^„(x) = f{x)e^^ + c.c, g|,^„(x) = ^f{x)e^7 + c.c. 

(which, with the notations of Theorem [T] corresponds to k = 1, 5 = and a = 
^i(k)a=(/o,i^/T)- 

The asymptotic models derived in Section [2] read in this particular case as follows: 

• The full dispersion model. The exact solution to (|4ip is approximated by 

<pp,i(^'^) = U(i)it,x)e'^^^-^ +C.C., 
where C/(i) — (/(i),(7(i)) solves (recall that D — —idx) 

(43) I + ^7T^Sf?wi^(i) " 7|l/(i)l'/(i), 
I 5(1) = "^/(i) (polarization condition). 

• The nonlinear Schrodinger equation. The exact solution to (j4ip is approx- 
imated by 

ulpp,2(t^ = U{2) (t, x)e''^-^ + c.c, 
where U(2) = (/(2),g(2)) solves 

1 + - |9^/(2)) = £^l/(2)P/(2), 

1 5(2) — ^^/(2) (polarization condition). 

• The nonlinear Schrodinger equation with improved dispersion relation. We 
approximate the solution of (|^T|) by 

Kpp^sit^x) = U(^3){t,x)e'^^ +C.C., 
where C/(3) = (/(3),g(3)) solves 



(44) 



(45) 



(1 - ted, - e'^dDdJ^,-, + 72(^- - + ^'^^')/(3) = 7|I/(3)IV(3), 
(7(3) = "^/(s) (polarization condition); 

note that this modified Schrodinger equation corresponds to the set of co- 
efficients (HOl). 
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3.2. The numerical scheme. We use a spectral method m space and a splitting 
technique in time for all the equations introduced in the previous section. We give 
here some details on the numerical scheme used for (|4ip : for ([33]), and pS)) . 
we use straightforward adaptations of this scheme. 

Let us denote by SL{t) and SNLit) the evolution operator associated respectively 
to the linear and nonlinear part of ()4ip : namely, 




SL{t)u° ^ if (t), git)), with { f and (/(0),.g(0)) 



and 

SNLit)u' = if it), git)), with I a;^r;(|)p+|^|f)P^ and (/(0),ff(0)) = u° 
(and with periodic boundary conditions). 

The numerical computation of 5*^ (t) is made through an FFT-based spectral method 
while an explicit integration is used for SNL{t)', we then use a second order split- 
ting scheme to compute m""*"^ ~ + l)At) in terms of u" ~ u(nAt) (where At 
denotes the time step): 

= SLi^)SNLiAt)SLi^)u^\ 

3.3. Numerical results for short pulses. In this section, we are interested in 
short pulses, that is we consider initial conditions for (|¥T|) of the form with 

/°(x) = G(^), 

where G is a smooth function. In the present numerical computations, the compu- 
tational domain is [0,L] with L = SOtt, and we take xq = 15 and G(a;) = e^^ . 

The accuracy of the approximations ([43|l . ([44|) and (|45l) is checked using the 
following quantity: 

Uf^, F f m l/ftr) -(/(,) (^,-)e^"^+c.c.)|oo 

(46) S(,)(£,/3)= sup — — , 

te[o,i] |Jl,t, vloo 

where j = 1 for the full dispersion model ([321), J = 2 for the usual Schrodinger 
approximation (|44p . and ^ = 3 for our new modified Schrodinger equation (|45[) . 
The exact solution and the difference between the exact solution with the approx- 
imation furnished by the FD, Schrodinger and improved Schrodinger models are 
plotted in Figure [4] for e = 0.01, /3 = 0.075 and at time T = 50 on the domain 
X € [0, 307r]. The following computations are also performed to test the accuracy 
of the approximate models: 

• Test 1: With (3—1 fixed, we let e vary from e = 0.001 to e = 0.1. This configu- 
ration corresponds to usual wave packets for which the three models should have a 
comparable accuracy of G(e) when e is small enough. One can indeed observe on 
Figure [S] that the errors £'Q)(e,/3) (j — 1,2,3) grow hnearly with e. One will also 
check that when e is too large (e ~ 5.10^^ for a rough precision of 20%), none of 
the models furnishes a good approximation. 

• Test 2: Here, we look at the same configuration as in Test 1 but with (3 — 0.1, 
that is, we investigate here short pulses. We can observe on Figure [5] that the 
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Figure 4. Short pulses: the exact solution, and the difference 
between the exact solution and the FD, Schrodingcr and improved 
Schrodinger models (from left to right and top to bottom) with 
£ = 0.01, /3 = 0.075 and T = 50. 
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Figure 5. The errors Ei^j){e,f3) for /3 = 1 and e e [0.001,0.1]; 
j — 1 corresponds to FD, j = 2 to Schrodinger and j = 3 to the 
improved Schrodinger. 



FD and improved Schrodinger models provide a good approximation, but that the 
usual Schrodinger approximation is completely inaccurate. 

• Test 3: Here, e — 0.01 is fixed and we let /3 vary from (3 — 0.01 (short pulses) 
to /3 = 1 (wave packets). It can be checked that the FD model furnishes a correct 
approximation for (3 > 0.03 and that for such values of /3, the improved Schrodinger 
approximation has the same precision. This is to be contrasted with the usual 
Schrodinger approximation which is completely inaccurate until f3 ~ 0.2. 



3.4. Numerical results for chirped pulses. In this section, we are interested in 
chirped pulses, that is we consider initial conditions for (|411) of the form (|42| . with 



fix) = G{x - xo) cos(i cos(^^-— ^)), 
P P 
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Figure 6. Short pulses: The errors E(^j){e,l3) for /3 = 0.1 and 
e e [0.001, 0.1]; j = 1 corresponds to FD, j = 2 to Schrodinger and 
j = 3 to the improved Schrodinger. 
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Figure 7. Short pulses: The errors £^(j)(e,/3) for e = 0.01 and 
/3 e [0.01, 1]; j = 1 corresponds to FD, j = 2 to Schrodinger and 
J = 3 to the improved Schrodinger. 

where G is a smooth function. In the present numerical computations, the compu- 
tational domain is [0, L] with L = 307r, and we take xq = 15 and G{x) — e^^ . 

The accuracy of the approximations ((43|) . (I44p and (|45)) is checked using the 
quantities E(^j^{e,f3) {j = 1,2,3) defined in ([^ . 

The exact solution and the diffenrence between the exact solution and the FD, 
Schrodinger and improved Schrodinger models are plotted in Figure[5]for e — 0.01, 
/3 = 0.3 and at time T = 1/e = 100. The following computations are also performed 
to test the accuracy of the approximate models: 

• Test 1: With [3 — 0.1 fixed, we let e vary from e = 0.001 to e = 0.1. We can 
observe on Figure [H] that the FD and the improved Schrodinger models are good 
approximations for e < 0.003. Above this value, the approximation is no longer 
pertinent. Furthermore, the classical Schrodinger model is inapropriate for this 
range of parameters. 

• Test 2: Here, e = 0.01 is fixed and we let /3 vary from /3 — 0.01 (chirped pulses) 
to /3 = 1 (wave packets). We observe on Figure [TO] that both FD and improved 
Schrodinger models become appropriate for (3 > 0.1 whereas the Schrodinger ap- 
proximation is acceptable for (3 > 0.4. 
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Figure 8. Chirped pulses: the exact solution, and the difference 
between the exact solution and the FD, Schrodinger and improved 
Schrodinger models (from left to right and top to bottom) with 
e = 0.01, /3 = 0.3 and T = 100. 
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Figure 9. Chirped pulses: The errors E^jj{e,f3) for /3 = 0.01 and 
e G [0.001, 0.1]; j — 1 corresponds to FD, j — 2 to Schrodinger and 
j = 3 to the improved Schrodinger. 
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Figure 10. Chirped pulses: The errors i?Q)(e,/3) for e = 0.01 and 
/? G [0.01, 1]; j = 1 corresponds to FD, j = 2 to Schrodinger and 
j = 3 to the improved Schrodinger. 
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